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Epicyclic Frequencies and Resonant Phenomena 
Near Black Holes: The Current Status 
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We briefly discuss the theory of epicyclic motion near rotating black holes and its relation 
to the origin of high frequency Quasi-Periodic Oscillations (QPOs) seen in many cases 
of accreting black hole binaries. We also point out some new frequency ratios predicted 
by the theory. 

In the framework of general relativity the successive theory of epicyclic motion 
around rotating black holes was first developed in Refs. 1,2 (see also a review paper 3 ). 
As is known, the rotating black holes are uniquely described by the famous Kerr 
solution 

ds 2 = -^ (dt - a sin 2 d<P f + E ( ^ + d0 2 ) + ^ [adt - (r 2 + a 2 ) d<j>] 2 , (1) 

2—i \ ZA / 2-1 

where the metric functions 

A = r 2 + a 2 - 2Mr , E = r 2 + a 2 cos 2 9 . 
The equation of motion for a test particle of mass fi and charge e has the form 
(fx* dx« dx? e dx" 

ds 2 al3 ds ds fi u ds ' [ ' 

From the symmetry properties of the metric (JlJ it follows that a cyclic motion of 
the test particle is possible in the equatorial plane 8 = ir/2 . For the cyclic motion 
we have z fi (s) — u°s {1, 0, 0, u^}, where luq is the angular velocity of the motion. 

To describe the epicyclic motion of the test particle we employ a perturbative 
expansion about the cyclic orbits, i.e. we assume that a; M (s) = z^(s)+^(s). Substi- 
tuting this into equation ([2|) and adding also the external force / p = / A ' ly u ly (u ) _1 , 
we arrive at the equation 
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where A/" M stands for the non-linear terms, the quantities 7^ and d a U^ must be 
taken on a cyclic orbit r = r*o, 9 = n/2 and 

e ,.1 



(4) 



Performing an appropriate integration in ^ and restricting ourselves to the first 
order terms in £ we obtain the set of equations 

eft 
rf 2 C 



,l t +tf? = -jjjji I f A(tt > A = 0,3 = *, (5) 



r=-£o / r -^i // A A , (6) 



g?£ 2 r /zu 

We see that in the absence of the external force the epicyclic motion is governed 
by the two decoupled oscillatory-type equations. In other words, it amounts to free 
radial and vertical oscillations with the associated frequencies 

BIT 1 A 1/2 fdU e \ 1/2 

The explicit expressions for these frequencies were first found in Ref. 1 in the most 
general case of black holes with an electric charge or with an external uniform 
magnetic field. Further, the epicyclic frequencies were used to explore the maser 
effect near black holes 2 as well as the effects of periastron precession and the Lense- 
Thirring drag in the Schwarzschild and Kerr fields threaded by a cosmic string. 4 
For a "pure" Kerr field the frequencies are 1 ' 2 

2 2 / 6M 3a 2 \ 2 2 ( 3o 2 \ . . 

(J r = lj II — ± 8auj K I , uj e = lo I 1 + — =p 4aujK I , (9) 



where ujq = ±wj<-/(1 ± clujk), ^>k = y M/r 3 is the Kepler frequency and the two 
signs refer to direct and retrograte motions. Later, in 1987 our expression for the 
frequency of radial oscillations was re-derived in Ref., 5 while the formula for the 
frequency of vertical oscillations re-appeared 6 only in 1990. 

In Ref. 1 we have also pointed out that when taking into account the non- linear 
terms on the right-hand-side of Eq. Q a coupling between the oscillations occurs. 
If the frequencies u r and uig are in the rational relation k r uj r — kgtve, where k r 
and kg are integers, the resonances must take place. We have plotted the positions 
of low-order (k =\ k r \ + \ kg |) resonances: k = 3, 4, 5, or (2 : 1, 3 : 1, 3 : 2). In 
recent developments, this theory has been extensively used to explain the frequency 
ratios and the origin of high frequency Quasi-Periodic Oscillations (QPOs) observed 
in many cases of accreting black hole or neutron star binaries (see papers 7,8 and 
references therein). In some cases, even our plots for the positions of the non- linear 
resonances were exactly reproduced (see for instance, Ref. 9 ). 
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Our model also admits the appearance of the forced resonances under a per- 
turbing force. Solving Eqs.(5)-(7) for Fourier transforms we obtain 

t A( ^ i A £ r, \ e f A ( U ,™) / 1n x 

£ (u,m) = - — 7i £ (^ m )-T^o— . ( 10 ) 

tr< . e / r (^,m)-^ 7 3 l / A (w,m) 

fl e f e (uj,m) 
£ [u,m) = — — 5-, w m =w-mw - (12) 

We see that at the frequencies lo — mujQ , lu = Cl^r — mujQ ± uj r , lu — = 
mu>o ± ll>6 the system exhibits the resonant behaviour. 2 The similar type of the 
resonant frequencies was discussed by Kato 10 within a model for kHz QPOs reso- 
nantly induced in warped disks of X-ray binaries. In the latest developments, we 
have performed a complete numerical analysis of the frequency ratios pertaining to 
the orbits where the frequency of the radial oscillations approaches its maximum 
value. Some results for direct orbits and for masses ~ 10M© are given in the Table. 



a/M 


Trnax / 


^o(Hz) 


^y(Hz) 


i/«(Hz) 


vejv r 


vo/v r 




0.00 


8.00 


141.42 


70.71 


141.42 


2.00 


2.00 


1.00 


0.10 


7.58 


152.60 


75.73 


151.17 


2.00 


2.02 


1.01 


0.20 


7.15 


165.68 


81.52 


162.37 


1.99 


2.03 


1.02 


0.30 


6.70 


181.26 


88.27 


175.45 


1.99 


2.05 


1.03 


0.40 


6.24 


199.97 


96.29 


190.74 


1.98 


2.08 


1.05 


0.50 


5.76 


223.29 


105.99 


209.23 


1.97 


2.11 


1.07 


0.60 


5.26 


252.90 


118.03 


231.78 


1.96 


2.14 


1.09 


0.70 


4.71 


292.71 


133.53 


260.59 


1.95 


2.19 


1.12 


0.80 


4.11 


349.83 


154.57 


298.90 


1.93 


2.26 


1.17 


0.90 


3.42 


442.93 


185.95 


353.95 


1.90 


2.38 


1.25 


0.99 


2.45 


662.11 


243.45 


447.52 


1.84 


2.72 


1.48 
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